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We study the complete one loop contribution to W^V, V = Z,'y, both in the Minimal 

Supersymmetric Standard Model (MSSM) and in the Two Higgs Doublet Model (2HDM). We 
evaluate the MSSM contributions and compare them with the 2HDM ones taking into account 
6 — > 57 constraint, vacuum stability and unitarity constraints in the case of 2HDM, as well as 
experimental constraints on the MSSM and 2HDM parameters. In the MSSM, we found that in the 
[ ■ intermediate range of tan/3 < 10 and for large At, the branching ratio of — > W^Z can be of 

the order 10""^ while the branching ratio of — » W^y is of the order 10~^. We also study the 
effects of the CP violating phases of Soft SUSY parameters and found that they can modify the 
. branching ratio by about one order of magnitude. However, in the 2HDM where the Higgs sector is 

' less constrained as compared to the MSSM higgs sector, one can reach branching ratio of the order 

QQ I 10"^ for both modes. 
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I. INTRODUCTION 

The Standard Model (SM) of electroweak interactions Q is very successful in explaining all experimental data 
available till now. The cornerstone of the SM, the electroweak symmetry breaking mechanism, still has to be es- 
tablished and the Higgs boson has to be discovered. The main goals of future colliders such as LHC and ILC is to 
study the scalar sector of the SM. Moreover, the problematic scalar sector of the SM can be enlarged and some simple 
extensions such as the Minimal Supersymmetric Standard Model (MSSM) and the Two Higgs Doublet Model (2HDM) 
0, 1^ are intensively studied. Both in the 2HDM and MSSM the electroweak symmetry breaking is generated by 2 
Higgs doublets fields $i and <I>2. After electroweak symmetry breaking we are left with 5 physical Higgs particles (2 
charged Higgs iJ^, 2 CP-even H^, hP and one CP-odd A^). The charged Higgs , because of its electrical charge, 
is noticeably different from the other SM or 2HDM/MSSM Higgs particles, its discovery would be a clear evidence of 
physics beyond the SM. 

The charged Higgs can be copiously produced both at hadrons and e^e~ colliders. In hadronic machines, the 
charged Higgs bosons can be produced in many channels: i) the production of ti pairs may offer a source of charged 
Higgs production. If kinematicaly allowed mfj± < mt, the top quark can decay to H~^b, competing with the SM 
decay t — > W^b. This mechanism can provide a larger production rate of charged Higgs and offers a much cleaner 
signature than that of direct production. 

a) single charged Higgs production via gb tH~ , gg tbH~ , qb q'bH~ 0|. in) single charged Higgs production 
in association with gauge boson via gg — > W^H^ or bb W^H^ and also single charged Higgs production 
in association with boson via qq,gg A^H^ 0. iv) pair production through qq annihilation Q or gluon 
fusion. 

At e+e^ colliders, the simplest way to get a charged Higgs is through pair production. Such studies have 
been already undertaken at tree-level j^] and one-loop orders and shown that e^e~ machines will offer a clean 
environment and in that sense a higher mass reach. We mention also that charged Higgs bosons pair production 
through laser back-scattered 77 collisions has been studied in the literature ,10] and found to be prominent to 
discover the charged Higgs boson. 

Experimentally, the null-searches from L3 collaborations at LEP-II derive the lower limit of about m^± > 80 GeV 
[ni, a hmit which applies to all models (2HDM or MSSM) in which BR(i7± TVr)+ BR(i7± cs)^l. DELPHI 
has also carried out search for — > A^W^^ ^ topologies in the context of 2HDM type I and derive the lower limit 
of about mjLf± > 76 GeV Recently, CDF Run II excluded a charged Higgs mass in the range 80 < < 160 

GeV This limit can apply for the MSSM with low tan/3. If the charged Higgs decay exclusively to fi/ the 

BR(t — *■ 77+6) is constrained to be less than 0.4 at 95%C.L. On the other hand if no assumption is made on charged 
Higgs decay, the BR(i H^b) is constrained to be less than 0.91 at 95%C.L. 

At the LHC, the detection of light charged Higgs boson with m^± < nit is straightforward from top production 
followed by the decay t — *■ bH'^^ . Such light charged Higgs (m//± < mt) can be detected also for any tan/3 in the tv 
decay which is indeed the dominant decay mode {T^ . However, for heavy charged Higgs masses m^± > rrit which 
decay predominantly to tb, the search is rather difficult due to larg e irreducible and reducible backgrounds associated 
with — > tb decay. However, it has been demonstrated in [l3| that the tb signature can lead to a visible 

signal at LHC provided that the charged Higgs mass below 600 GeV and tan/? is either below < 1.5 or above > 40. 
Ref. proposed — > tv as an alternative decay mode to detect a heavy charged Higgs, even if such decay is 
suppressed for heavy charged Higgs it has the advantage being more clean than tb. 

An other alternative discovery channel for heavy charged Higgs is its decay to charged gauge boson and lightest CP- 
even Higgs: VF^/i", followed by the dominant decay of to bb [23|. Since the branching ratio of — > W^hf^ 
is suppressed for High tan (3, this channel could lead to charged Higgs discovery only for low tan/3 where the branching 
ratio of W^h^ is sizeable. 

Both in 2HDM as well as in MSSM, at tree level, the coupling — > W^j is absent because of electromagnetic 
gauge invariance J7(l)om- While the absence of —> W^Z is due to the isospin symmetry of the kinetic Lagrangian 
of the Higgs fields [21|. Therefore, decays modes like — > W^^, W^Z are mediated at one loop level and 

then are expected to be loop suppressed [2^ |2^ |2^ . We emphasize here that it is possible to construct models 
with an even larger scalar sector than 2 Higgs doublets, one of the most popular being the Higgs Triplet Model (HTM) 
(23. A noteworthy difference between 2HDM and HTM is that the HTM contains a tree level ZW^H^ couphng. 
Motivated by the fact that there is no detailed study about W^V, V = Z,j, in the the framework of MSSM 



^ Note that i n the 2HDM it may be possible that the decay channel —> W^A'^ is open and even dominate over tu mode for m£f± <, mt 
^ Note that at Tevatron run II, the charged Higgs is also searched in top decay Il6l . 
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in the literature which take into account b 57 and other electroweak and experimental constraints. We would like 
to reconsider and update the existing works [23, [2^ 123. 125I l26l | on the charged Higgs boson decays into a pair of 
gauge boson: i?^ ^ W^^i.W^Z both in 2HDM and MSSM with and without CP violating phases. Although these 
decays are rare processes, loop or/and threshold effects can give a substantial effect. Moreover, once worked out, 
any experimental deviation from the results within such a model should bring some fruitful information on the new 
physics and allow to distinguish between models. We would like to mention also that, those channels have a very 
clear signature and might emerge easily at future colliders. For instance, if W^Z is enhanced enough, this 

decay may lead to three leptons final state if both W and Z decay leptonically and that would be the corresponding 
golden mode for charged Higgs boson. 

Charged Higgs decays: W^j,W^Z, have received much more attention in the literature. W^Z 

has been studied first in the MSSM in Ref. [IJ has considered both W^-j and W^Z in the 

MSSM and show that the rate of W^j is very small while the rate of W^Z can be enhanced by heavy 

fermions particles in the loops. The four th g eneration contribution was given as an example. Although the squarks 
contribution has been considered in Ref. |23|, Left-Right mixing which could give substantial enhancement has been 
neglected. — *■ W^^ was also studied in [23| within the MSSM. Later on, Ref. studied the possibility of 
enhancing ~> W^Z by the non-decoupling effect of the heavy Higgs bosons in the context of 2HDM, substantial 
enhancement was found j^^. Recently. — > W^j was also studied in 2HDM t ype H [2^. All the above studies 
has been carried out either in unitary gauge [22, 1131 or in the nonlinear i?j-gauge |26l | . The analysis of and |23 
have been performed in 'tHooft-Feynman gauge without any renormalization scheme. It has been checked in l25l | 
that the sum of all Fcynman diagrams: vertex, tadpoles and vector boson-scalar mixing turns out to be Ultra- Violet 
finite. 

In the present study, we will still use 'tHooft-Fcynman gauge to do the computation. However, the amplitudes of 
ij± ^ w^j and ^ W^Z are absent at the tree level, complications like tadpoles contributions and vector 
boson-scalar mixing req uire a careful treatment of renormalization. We adopt hereafter the on-shell renormalization 
scheme developed in j28| . 

The paper is organized as follows. In section II, we describe our calculations and the one-loop renormalization scheme 
we will use for W^Z and — > W^^. In Section III, we present our numerical results and discussions, and 

section VI contains our conclusions. 



II. CHARGED HIGGS DECAY: H± ^ VV^V 



As we have seen in the previous section, both in 2HDM as well as in MSSM, at tree level, the coupling 
W^^ and ff^ — *■ W^Z do not exist. They are generated at one loop level and then are expected to be loop 
suppressed [l^, 113, 0, HE US • Hereafter, we will give the general structure of such one loop couplings and discuss 
the renormalization scheme introduced to deal with tadpoles and vector boson scalar boson mixing. 



One loop amplitude W^V 



The amplitude M for a scalar decaying to two gauges bosons Vi and V2 can be written as 



3 Li* 7/ * 



where ey^ are the polarization vectors of the Vi. 

According to Lorcnz invariance, the general structure of the one loop amplitude A^^^^ of S* — ^ V^VJ decay, if CP is 
conserved, is 

(2) 

where pi^2 are the momentum of Vi, V2 vector bosons, .^1^2, 3 are form factors, and e^vpa is the totally antisymmetric 
tensor. The form factor J-i has dimension 2 while the other are dimensionless. The analytic expression for J-i are 
given in appendix B. 

For — !■ W^j, electromagnetic gauge invariance implies that =l/2{?n'^ — ■m'jj^)J^2 ^3- This means that 
only T2 and will contribute to the decay W^-f. In case of W^Z, there is no such constraint on 

form factors. 

In terms of an effective Lagrangian analysis, from gauge invariance requirement we can write: 

Ceff = giH^W;^V'' + 52i^*i^^^M/^.. + igse^.paH^ F^" F'^p^ + h.c (3) 
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the first operator H^W^V^ is dimension three and the last two operators H^Fy'Fwuv and e^^p^j jj^p^^^pWpa 
are dimension five. One conchide that 52, 3 (resp. 171) must be of the form g{R)/M (resp Mg{R)) with M a heavy 
scale in MSSM or in 2HDM. g{R) a dimensionless function and i? is a ratio of some internal masses of the model 
under studies. Therefore, it is expected that in case of ^W^Z decay, will grow quadratically with internal 
masses while J^2.3 will have only logarithmic dependence. A contrario, for — > W^j decay, the electromagnetic 
gauge invariance relates JFi and J^2 and then the amplitude of — > W^'f will not grows quadratically with internal 
masses. One expect that the decay — > W^j is less enhanced compared to W^Z. 

After squaring the amplitude and summing over polarization vectors, the decay widths as functions of the form factors 
J- 1 , Ti and take the following form: 
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where A = (m^± - m^^ - ,,^2; -■t,,,^,,,^ 



B. On-shell renormalization 



We have evaluated the one- loop induced process — > W^V in the 'tHooft-Feynman gauge using dimensional 
regularization. The typical Feynman diagrams that contribute to — > W^V are depicted in Fig. 1. Those diagrams 
contains vertex diagrams (Fig. 1.1 1.11), W^-H^ mixing (Fig. 1.12 1.14), H^-G^ mixing (Fig. 1.15 1.17) 
and H^-W^ mixing (Fig. 1.18 1.20). 

Note that the mixing H^-W^ (Fig .1.12, 1.13, 1.14) vanishes for an on-shell transverse W gauge boson. There is 
no contribution from the W^-G^ mixing because •yG^H^ and ZG^H^ vertices are absent at the tree level. All the 
Feynman diagrams have been generated and computed using FcynArts and FormCalc 12^ packages. We also used 
the fortran FF-package [s^l in the numerical analysis. 

Although the amplitude for our process is absent at the tree level, complications like tadpole contributions and vector 
boson-scalar mixing require a careful treatment of renormalization. We adopt, hereafter, the on-shell renormalization 
scheme of [sJl, for the Higgs sector, which is an extension of the on-shell scheme in [s^- In this scheme, field 
renormalization is performed in the manifest-symmetric version of the Lagrangian. A field renormalization constant 
2 is assigned to each Higgs doublet $1,2 • Following the same approach adopted in the Higgs fields and 
vacuum expectation values Vi are renormalized as follows: 



(6) 



With these substitutions in the scalar covariant derivative Lagrangian of the Higgs fields (in the convention of |3|), 
followed by expanding the renormalization constants Zi = l+SZi to the one-loop order, we obtain all the counter-terms 
relevant for our process: 



6[Z,W^H^] 

where k denotes the momentum of the incoming and 



Sw 



Cw 



sin 2/3 5v2 5vi 
A = — - — mw[ 

2 V2 Vl 



SZip^ — 6Z^^] . 



(7) 
(8) 

(9) 



(10) 



Denoting the one particle irreducible (IPI) two point function for W^H^ (resp G^H^) mixing by ±ifc^I]^^±jLf± (fc^) 
(resp iT,Q±H± (k^)) where k is the momentum of the incoming (resp G*), and is outgoing. The renormalized 
mixing will be denoted by E. 

In the on-shell scheme, we will use the following renormalization conditions: 
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FIG. 1: Generic contributions to W^V 



• The renormalized tadpoles, i.e. the sum of tadpole diagrams u and tadpole counter-terms ?>h,H vanish: 

Th + &th = 0, Th + = . 

These conditions guarantee that v\^i appearing in the renormalized Lagrangian are located at the minimum 
of the one-loop potential. 

• The real part of the renormalized non-diagonal self-energy S//±vi/±(fc^) vanishes for an on-shell charged Higgs 
boson: 

5ReS^±vi/±(m2^±) = (11) 
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FIG. 2: Generic contributions to — > and mixing as well as counter-terms needed. 



This renormalization condition determines the term A to be 
and consequently 5[A^W^H^] and (5[Z^W±i/=F] 

are also fixed. 



(12) 



The last renormalization condition is sufficient to discard the real part of the H^-G^ mixing contribution as well. 
Indeed, using the Slavnov-Taylor identity [s^ 



fc^E^±i4/± (fc^) — TOvyE//±G± (fc'^) = at = m\j± 
which is valid also for the renormalized quantities together with eq. it follows that 

afJeS/fi G± (m^± ) = . 



(13) 



(14) 



In particular, the Feynman diagrams depicted in Fig. 1.9 will not contribute with the above renormalization conditions, 
being purely real valued. 

To make the amplitude of Fig.l Ultra- Violet finite we need to add the following counter-terms: counter-terms for 
^W^H^ and ZW^H^ vertices Fig. 2. 2a, a counter-term for the W^-H^ mixing Fig. 2. 2b, 2. 2d, and a counter-term 
for the G^-H^ mixing Fig. 2.2c. 



III. NUMERICS AND DISCUSSIONS 



In our numerical evaluations, we use the following experimental input quantities |3J|: a ^ = 129, mz, mw, mt, 
nib = 91.1875, 80.45, 174.3, 4.7 GeV. In the MSSM, we specify the free parameters that will be used as follow: i) The 
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FIG. 3: Branching ratios of W^Z (left) and iJ± -> W±7 (right) as a function of m„± in the MSSM and 2HDM for 

MsusY = 500 GeV, M2 = 175 GeV, /i = -1.4 GeV and At = At = Ar = -^t for various values of tan/3. 
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FIG. 4: Branching ratios for H 
GeV, M2 = 200 GeV, m^± = 500 GeV, At 



W^Z (left), 



At 



W^-y (right) as a function of tan /3 in the MSSM with Msusy = 500 
Ar = — M for various values of At- 



MSSM Higgs sector is parametrized by the CP-odd mass rriAa and tan /?, taking into account radiative corrections 
from |37l |. and we assume tan/3 > 3. ii) The chargino-neutralino sector can be parametrized by the gaugino-mass 
terms Mi, M2, and the Higgsino-mass term fjL. For simplification GUT relation Mi k, M2/2 is assumed, iii) Sfermions 

are characterized by a common soft-breaking sfermion mass Msusy = ^J^l = Mji, fi the parameter and the soft 
trilinear couplings for third generation scalar fermions At^t.r- For simplicity, we will take At = = A^. 

When varying the MSSM parameters, we take into account also the following constraints: i) The extra contributions 
to the Sp parameter from the Higgs scalars should not exceed the current limits from precision measurements [s^ : 
\Sp\ < 0.003. ii) b ^ constraint. The present world average for inclusive b ^ rate is [s^l 13{B Xsj) = 
(3.3 ± 0.4) X 10"''. We keep the B Xsj branching ratio in the 3cr range of (2.1-4.5) x 10"''. The SM part of 
B Xsj is calculated up to NLO using the expression given in While for the MSSM part, the Wilson coefficient 
C7 and Cs are included at LO in the framework of MSSM with CKM as the only source of flavor violation and are 
taken from js^. iii) We will assume that all SUSY particles Sfermions and charginos are heavier than about 100 
GeV; for the light CP even Higgs we assume m/jO > 98 GeV and tan/3 > 3. 

The total width of the charged Higgs is computed at tree level from Q without any QCD improvement for its 
fermionic decays iJ± //'. The SUSY channels like H+ f,f^ and H+ — !■ XiXj' ^I'G included when kinematicaly 
allowed. In Fig.|2l we show branching ratio of — > W^Z (left) and W^j (right) as a function of charged 

Higgs mass for tan/3 = 16 and 25. In those plots, we have shown both the pure 2HDM (only SM fermion, gauge 
bosons and Higgs bosons with MSSM sum rules for the Higgs sector) and the MSSM (2HDM and SUSY particles) 
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FIG. 5: Scatter plot for branching ratios of W^Z (left), 77='= -+ W^y (right) in the {■mH±, tan/3) plane in the MSSM 

for MsusY = 1 TeV, Ma = 175 GeV, At,b,r = Msusy and = -1 TeV. 



10^ 

N 

+1 

+1 

I 
m 



10' 



tanS = 3 



..IQ... 



= 20 



10"^ 



200 400 600 800 1000 1200 1400 1600 
A,(GeV) 



Tio" 

+1 

X 



10"' 



tanP = 3 


1 






\ k 





200 400 600 800 1000 1200 1400 1600 
A,(GeV) 



FIG. 6: Branching ratios for ^W^Z (left) and W^-y (right) as a function of At in the MSSM with Msusy = 500 
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contribution. As it can be seen from those plots, both for W^Z and W^j the 2HDM contribution is 

rather small. Once we include the SUSY particles, we can see that the Branching fraction get enhanced and can reach 
10~ in case of W^Z and 10"^ m case of W^j. The source of this enhancement is mainly due to the 

presence of scalar fermion contribution in the loop which are amplified by threshold effects from the opening of the 
decay — > tib*. It turns out that the contribution of charginos neutralinos loops does not enhance the Branching 
fraction significantly as compared to scalar fermions loops. The plots also show that, the branching fraction is more 
important for intermediate tan/3 = 16 and is slightly reduced for larger tan/3 = 25. 

This tan/3 dependence is shown in Fig. ^ both for W^Z and — > W^j for three representative values of 

At- It is obvious that the smallest is tan/3 the largest is the branching fraction. Increasing tan/3 from 5 to about 40 
can reduce the branching fraction by about one or two order of magnitude. 

We also show a scatter plot Fig. O for W^Z (left) and W^j (right) in (m//±, tan/3) plane for 

At = — /i = 1 TeV, Msusy = At and M2 = 175 GeV. As it can be seen from Fig [31 there is only a small area for 
tan/3 < 10 where the branching ratio of W^Z can be in the range 10~^-10^'^. 

We now illustrate in Fig. El the branching fraction of — > W^Z (left) and (right) as a function of 

At = At = Ar — — /i for Msusy = 500 GeV and M2 — 200 GeV. Since b ^ favor At and /i to have opposite 
sign, wc iix fi ~ —At and in this sense also fi is varied when At is varied. Both for — > W^Z and W^j, 
the chargino-neutralino contribution which is rather small decrease with fj, = —At, the largest is At the smallest is 
chargino-ncutralino contribution. In case of W^Z, for At < 1 TeV it is the pure 2IIDM contribution which 
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dominate and that is why it is almost independent of At while for large At the branching ratio increase with At- It is 
clear that the largest is At the largest is the branching ratio which can be of the order of 10^"^ for W^Z with 

tan /3 — 10. As wc know from h° ~^Sn ^"^^ ^ MSSM [s^, the squarks contributions decouple except in the 

light stop mass and large At limit jS^. In W^V case, the same stuation happen. As we can see from Fig.El 

(left), for intermediate At, 300 < At < 1000 GeV, the squarks are rather heavy and hence their contributions is small 
compared to 2IIDM one. While for large At the stop becomes very light < 200 GeV and hence enhance W^V 
width. Of course this enhancement is also amplified by H^tL,Fib*ji and H^tl^jiVtl couplings which are directly 
proportional to At^b,T- In case of W^j decay, the pure 2IIDM and sfermions contributon are of comparable 

size, the branching ratio increases with At- 

We have also studied the effect of the MSSM CP violating phases. It is well known that the presence of large 
SUSY CP violating phases can give contributions to electric dipole moments of the electron and neutron (EDM) 
which exceed the experimental upper bounds. In a variety of SUSY models such phases turn out to be severely 
constrained by such constraints i.e. Arg(/x) < (10~^) for a SUSY mass scale of the order of few hundred GeV pOl| . 
For — > W^Z and — > W^^ decays which are sensitive to MSSM CP violating phases through squarks and 
charginos-neutralinos contributions, it turns out that the effect of MSSM CP violating phases is important and can 
enhance the rate by about one order of magnitude. For illustration we show in Fig. [7| the effect of At^b,T CP violating 
phases for Msusy = 500 GeV, At^b,T = —j-i- = 1 TeV and M2 = 150 GeV. For simplicity, we assume that /i is 
real. As it is clear, the CP phase of At^b,T can enhance the rates of both W^Z,j by more than an ordre of 

magnitude. Those CP violating phases can lead to CP-violating rate asymmetry of decays, those issues are going 
to be addressed in an incoming paper [4ll| . 

We now turn to discuss the pure 2IIDM contribution to — > W^Z,W^^. For the 2IIDM parameterization we 
follow closely the notation of |43|. In our discussion we will take as free parameters: 

A5 , m/,0 , rriHO , to^o , mH± , tan/? and sin a (15) 

a is the CP-even mixing angle and A5 is the term that breaks softly the discrete symmetry $i — (f>i in the 2IIDM 
Lagrangian. 

To constrain the scalar sector parameters we will use both vacuum stability conditions as well as tree level unitarity 
and 6p constraints. In our study, we use the vacuum stability conditions from |43l | and tree level unitarity constraints 
are taken from 01 • 

As B ^ Xs"f constraint is concerned, it has been shown in that for models of the type 2IIDM-II, data on B ^ Xs"f 
imposes a lower limit of m//± > 350 GeV. In type I 2IIDM, there is no such constraint on the charged Higgs mass 
|38j . In our numerical analysis we will ignore these constraints and allow mH± Si 200 GeV in order to localize regions 
in parameter space where the branching ratios are sizeable. 

In the 2IIDM, the source of enhancement in — > W^Z and W^j decays can be either from bottom 

Yukawa coupling which is enhanced at large tan/3 or from the Higgs sector contribution through diagram like Fig. 1.2 
with (S',,S'j,5'fc) = (i?±,/i'',7J±) or (iJ±,ffO,ir±) and Fig. 1.8 with {S,, Sj)={H^ , h°) or {H^,H°). In contrast to 
the MSSM where the trilinear scalar couplings h^H^H~ and H^H~^H~ are function of the gauge couplings only, in 
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FIG. 8: The branching ratios of W^Z (left) and VV^-y (right) in type-II 2HDM model in (m^±,tan/3) plane for 

100 GeV < nii,o < 130 GeV, 150 GeV < hiho < 350 GeV, 100 GeV < mj^o,n± < 500 GeV, -1 < sina < 1, 0.1 < tan/? < 10, 
and As = 0. 
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FIG. 9: Branching ratios ff± W'^Z as the function of the tan (3 in 2HDM-II for m^o = 100 GeV, m^o = 200 GeV, 
m^o = 300 GeV a = P — n /2 and various values of A5. Left plot is for m^fi = 180 GeV and right one is for m^fi — 220 GeV 

the 2HDM those couplings are function of Higgs masses, a, tan/3 as well as A5 as it can be seen from their analytic 
expressions (in the notations of Rcf. ) 

H^H^H^ = ^— — - (m?ro (cos^ /3 sin a + sin^/3cosa) + m?r± sin2/3cos(/3 - a) - sin(/3 + a)\^v^) 

mw sm 2p 

h°H^H~ = — (m?o(cosacos^/3 - sinasin^/S) + m?,± sin2/3sin(/3 - a) - cos(/3 + a)\v^) (16) 

mw sin 2/3 

Even after imposing unitarity and vacuum stability constraints, those couplings can get large values compared to their 
MSSM values and that is the main difference between 2HDM scalar sector and the MSSM one. 

We have first compared our results to previous one given in Re fs. | 25l l26l| and found perfect agreement. However 
in Refs. |25ll26j . unitarity constraints were not imposed while Ref. |26| did not consider the case of W^Z. We 

have performed a systematic scan over the 2HDM parameters and found that both W^Z and W^^ 

can reach 10~^ branching ratio. In FiglHlwe have illustrated our results as scatter plot in the (tan/3, m//±) plane. As 
one can see from the plot, just before the opening of tb decay, for m^± < rrtf + mi,, the branching ratio of 

— !■ W^Z and — > W^^ can be of the order 10~^-10^'^. Once the decay H~ tb is open for mjj± ^ + mt, 
the charged Higgs width becomes large and the branching ratio of — > W^Z and W^j are suppressed. 

However as one can see, there is still a large area in the (tan /3, m[j± ) plane where the branching ratio of W^Z 
and W^^ can be larger than 10"'*'. 
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We illustrate in Fig. the branching ratio of W^Z as a function of tan/3 for various values of A5 and for 

= 180 GeV (left plot) and mfj± = 220 GeV (right plot). It is clear from both plots that the branching ratio is 
slightly enhanced for vanishing A5 case, this effect has been noticed also by In the left plot we are very close to 
H~ tb threshold, the branching ratio is larger than 10""^. Away from the threshold H~ — > tb and for m^± = 220 
GeV (right plot) the branching ratio of if^ — > W^Z is slightly reduced to the level of 10^**. 



IV. CONCLUSION 



In the framework of MSSM and 2IIDM we have studied charged Higgs decays into a pair of gauge bosons namely: 
if± W^Z and W^^j. In the MSSM we have also studied the effects of MSSM CP violating phases. In 

contrast to previous studies, we have performed the calculation in the 'tHooft-Feynman gauge and used a renormal- 
ization prescription to deal with tadpoles, W^-H^ and G^-H^ mixing. The study has been carried out taking 
into account the experimental constraint on the p parameter, b sj constraint, unitarity constraints, and vacuum 
stability conditions on all scalar quartic couplings Xi in 2IIDM case. Numerical results for the branching ratios have 
been presented. In the MSSM, we have shown that the branching ratio of — *■ W^Z can reach 10"'^ in some cases 
while 7 never exceed 10 ^. The effect of MSSM CP violating phases is also found to be important. 

In the 2IIDM we emphasize the effect coming from the pure trilinear scalar couplings such as hPH^H^ and H^H^ H~ . 
We have shown that in 2HDM both -*W^Z and 7 can have a branching ratio in the range 10 ^-10 ^. 

Those Branching ratio in the range 10^^-10"^ might provide an opportunity to search for a charged Higgs boson at 
the LHC through W^Z. 
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APPENDIX A: LAGRANGIAN AND COUPLINGS 



In this appendix, we list our notation for the gauge and Yukawa couplings used in this section. 
The photon coupling to fermions/sfcrmions are: 

^7// = -^f^w, g^lJ^ = -ejswSij (Al) 

The Z boson coupling to fermions/sfcrmions are: 

9zff = i-T3 + efs'^)/cw, g§ff = efs'^/cw (A2) 
9znf, = c^j/cw (A3) 



where the couplings Cij are given by 



W 



-2^3 sin 26*/ T3 sing^ -e/s 



(A4) 



where T3 = 1/2 for up quarks and T3 ~ ^1/2 for down quarks and Of is the mixing angle in the sfermion sector 
defined by: 

where /l.r are the weak cigenstatcs and /1.2 arc the mass eigcnstates. 
The W gauge boson coupling to a pair of fermions/sfcrmions is: 

9wff' = ' 9^ff' - (A6) 

9wfJ, = -^^(Xi (A7) 



The chargino mass matrix is: 



X 



M2 \f2mw sin /? 
\/2mw cos fi /i 



which is diagonahzed by the unitary matrices U and V via VX'^'U'^ = Md. 
The neutraUno mass matrix is: 



Y = 



/ Ml —mzsw cos f3 mzsw sin f3 

M2 mzcw cos (3 —mzcw sin (3 

—mzs^v cos P mzCw cos (i — /x 

\ mzsw sin P —mzCw sin (i —fi 



which is diagonahzed by the matrix TV via NY'^N'^ = No. 
The matrices that enter the W^XiXj couphngs are defined as: 



Of 



The matrices that enter the Zxfxf couphngs are defined as: 



iL 



2 



1. 



2 

W 



The matrices that enter the ZxiXj couphngs arc defined as: 

1 



The matrices that enter the H^x^xf ^^^d H^Ujdj couphngs are defined as: 



tL 



N*^v;^ + -j={n:^ + N*, t^ndw)v;^, 



1 



L,R 
3H±Uidi 



{nidi tan /3, rnuj / tan (3) ^ n _ , nid, 

' 3G±Uidi 



raw 



^/2l 



mw 



The couphngs W^G^V, W^W^V and H^H^V are given: 





5wGy 


gwwv 


9HHV 


z 


-mzs'^ 


Cw 






7 


mwsw 


-sw 




Sw 



The couphng ff are given by 

dH-bLtL 

dn-bLtR 

dH-bntL 



—mw sin 2P/V2 
—mw sin 

—mw sin 2/3/ V2 + (ml tan f3 + m^ cot /3) /V2mw 
mt{ii + At cot /3)/-\/2mn^ 
mf,{fJ. + Abtan (3) / V2mw 
mtmf,(tan p + cot P)/V2mw 



dHbiii 
9Hbii2 

9Hb2i2 



setSe.gH-bjL ~ '^StSe.gH-br.in - ^etCe,gH-bRh +C0tCe^gH-b 
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APPENDIX B: MSSM CONTRIBUTION FOR W^V 

We now list our results for the MSSM diagrams. The matrix elements for W^V {V = Z, 7) in the MSSM are 

partially presented here. We give only the fermionic contribution (SM fermion and chargino-neutralino contribution) 
as well as the scalar fermion contributions. For the bosonic contribution (higgs bosons and gauge bosons) we refer to 

.....+ 

It is to be understood that diagrams involving charginos X; f'i'o summed over i = 1,2 and diagrams involving 

neutralinos Xi ^re summed over i = 1, . . . , 4. 

The fermion triangle that enters Fig. ^la and Fig. ^Ib was computed in Ref. and agrees with our result. 

Fig. ^Ic with Xi, X? s-nd xt in the loop is analogous to the top/bottom quark triangle diagram and can be checked 
by substituting top /bottom quark couplings into the gaugino couplings. 
Diagram Fig. la: In the convention of cq[21 we have 



Ti = Ncgw{'2'{mt9vbb9H + m.bgvbb9H - mbgvbb9H)Bo - 2mtmbgvbb{^tgH + mbgH)Ca 



(Bl) 



- mbgybbgni'm-Hi: + f^w - ™y) " mtgvbb9H{--m, 



H± 



3m 



w 



'^'mbgvbb9H'm^ 



T2 



mg^bb9^{Smj,, + - ml) + m,.g^,,.g^«± + - m^) - 2g^,,gj,m],,m, 
■^gvbbinT-tg^ + 'mbgH)Coo^ 

-2Ncgw{jn-tgvbb9H Co - {m-tgybbgH - mbgybbgn) Ci + 3gyf,bimtgH + mbgn) C2 
"^NcgwijntgvbbgH Co + {mtgvbb9H + mbgy^hgH) Ci + gybbi'^tgH + mbgn) C 



(B2) 
(B3) 



Where Nc = 3 for quarks and 1 for leptons. The arguments of the Passarino-Veltman functions Bi and Ci are 
Bi{ml, ml, ml), C.i{m^,m1, m\j± , ml, ml, ml). 



Diagram Fig. lb: 

Similar to Fig.ni(a) with the exchange of 



mt ^ mt, 



R L L L 

gvbb, gvtt ^ gvbb, gn 



gn, gw^ gw 



Diagram Fig.QJlc 



(B4) 



^1 = 



{gmgb + 523455) + "ixj (515355 + 525455) - (5154.95 + 525355) 

™Xi(~2.9i54g5TOxi»7ixfc - "^gmgenixjmx^^ + 2gig3gemx,mx^ + 2gig:ig^m\ 



Bn 



+51.9355 (™^± + -i^v)) 



Cn 



- (515455 + 525355) ("i^± + TO^ - ml) 



+'™x. (515355 - 525455) (-?Ti|f± - 3m^ + my) - 2(3133,96 - 525455)"^xj"Ih' 



(515355 + 525456)(3to^± + - my) + m^^ (525455 + 5i5355)("t^± + - my) 



-^gmgs - 525355 )TOxfc'^ff± 



"^x. (515355 + 525455) + (515355 + 525455) 



c, 



00 



^2 



^3 = 



2( ("^xi (515355 + 525455)) Co - m^^^ (52.9356 - 51.9455) + "i^, (5i5355 + 525455) 



Ci 



Cs 



"^x. (515355 + 525455) + m^j (515355 + 525455) 
8(nT-Xi (515355 - .925455) Co + m^^ (525355 - 5i5455) + "^x. (5i5355 - 525456) 
m-xi (515355 - 525455) - mxj (52.9455 - 515355) C2) 



Ci 



(B5) 



(B6) 



(B7) 



The arguments of and Ci functions arc Bi{m1,m^.,m'^.), Ci{m'y^r,my,m^^±,m'^.,m'^^,m'^^). The couplings are 
given in the following table: 
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Xi Xj Xk 


9i 


92 


93 


54 


55 


96 


x° x° x+ 

^+ ^0 


Of, lew 


-^O'l^/cw 




^ki 







Diagram Fig. ^2: 



^2 = Nc gXg^gf (Co + 2(Ci + 2Ci2 + C2)) 



(B8) 
(B9) 



Where Nc ~ 3 for scalar quarks and Nc = 1 for scalar Icptons. the arguments of Ci and Cy functions are as 
Ci,ij{rn'^ , rriy , ni^,, rn^^, rri^,) . The couplings are given in the following table: 







Sk 


7 

5i 


g? 


g? 


w 

53 




vi 


II 


g-fvi' 


—gzvv 


gHui 


gwvi 


II 


II 




g^ii 


gzu 


gmv 


gwiv 


UL 


UL 


dL 


g^uu 


—gzuu 


gUud 


gWiid 


dL 

u 
k 


dL 

ij 
h 


UL 

bk 

ik 


g^dd 

r,uL ryt ryt ^uR ryt ryt 
57 ^li^lj g-y ^2i^2j 
^dL r)b TDh 1 ^dR r)b Tib 

g-i ^li^ij + g-r ^2i^2j 


gZdd 

r,uL -Dt rit n^R l?t 

gz ^li-^lj gz ^2i^2j 
^dL TDh T)b 1 ^dR d6 d6 
5z Rli^lj + gz R2iR2] 


gUdu 

gmi 
gmi 


gwdu 

r>t jyb 
^Ij^lk 
Tit r)b 

^Ik-^lj 



where = = cos9f and —R12 = R21 = sinOf. 
Diagram Fig. ^8: 

Ti = -Ncgig2Bo{m%±,ml.,ml.) 
The couplings are given in the following table: 



Si Sj 


gi 


92 




9H±uih 


guhw^v 


h II 


9H±iLbL 


giLbLW±V 


ii h 


9H±iibi 


5ti&iW±V 


ii 62 


9H±iib2 


9iJ>2W±v 




5if±t2bi 


9t2biW^V 




9H±iib-, 


guboW^V 



Diagram Fig. [212: 
:F2 = - 



2^c53 



mp, (5154 + 5255)^0 + (5154 + g2gb)mFi + (5254 + gigb)mFj Si) 



The arguments for the Passarino-Veltman functions Bi are Bi(m'^,rn'^p,,rn'^p, 
The couplings are given in the following table: 



Fj Sk 


51 


52 


53 


54 


55 


t b 

X° X+ 





9w±tb 


gH^w^v 
gH±w±v 


9H±tb 


9Hftb 



Diagram Fig. ^ 13: 



T2 



Nc 51 52 .93 



The couplings are given in the following table: 



(BIO) 



(Bll) 



(B12) 
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c 




Q 

Ok 


91 


92 


93 


Vl 


1 


n± 
ri 




9H±i>ilL 


9i>iiLW± 


ul 


d-L 


u± 


9h±h^V 


9H±U[^dL 


9uLdLW± 


ii 


bi 


ij± 


9H±H=fV 


9H±iibi 


9iibiW± 


ii 




H± 


9h±h=fv 


9H±iib2 


9iib2W± 




h 




9h±h^v 


9H±i2bi 


9i2biW± 




h 




9h±h^v 


dH^i^b-, 


9i.>b,W± 



Diagram Fig. 2)15: 



^1 



, (.9253 + .9154) Ao{m%^ ) + (5233 + 9i9i)mF, + {gi93 + 9294)mF^ 

+ [9293 + gigA)gbm]j±Bi 

The arguments of Bi functions are Bi , m|n. , m^p.) . The coupHngs are given in the following table: 



(B13) 



Fi Fj Sk 


51 


92 


93 


54 


55 


t b G± 
X° X+ G± 


9G±tb 


9G±tb 


9g±w±v 
9g±w^v 


9H±tb 


9HTtb 



Diagram Fig. ^16: 



^ ^C5i.92.93 „ , 2 2 2 N 
mj^^ — rrig ^ 



(B14) 



The couplings are given in the following table: 



5. 


5, 


Sk 


51 


92 


53 




h 


G± 


9H±i>ilL 


5g±w±v 


9G±i>ih 


UL 


dL 


G± 


9H±ULdL 


50*1^*^ 


9G±ULdL 


ii 


bi 


G± 


9H±iibi 


5g±w±v 


9G±iibi 


ii 


62 


G± 


9H±iib2 


5g±w±v 


5G±t 162 


t2 


h 


G± 


9H±i2bi 


5G±H'±y 


9G±i2bi 


t2 


b2 


G± 


9H±Ub,. 


5g±w±v 


9G±Ub2 



Diagram Fig. 18: 
^2 = 
^1 = 



Nags 



"Hi: - 



(^0(5 



254 + gigbjmFi 



(5254 + gigh)mF, + (5154 + g2gb)mF^ 



(B15) 
(B16) 



The arguments of Bi functions are as Bi {m^^ , m^, , m?p. ) . The couplings are given in the following table: 



F^ F, Vk 


51 


52 


53 


54 


55 


t b 


9w±tb 





9w±w±v 
gw±w±v 


9H±tb 


„i?, 

9H±tb 



Diagram Fig. ^ 19: 



T2 



{my - m\f)T2 
515253 



.2 ^2 



Bo + 2Bi 



(B17) 
(B18) 



The arguments of Bi functions are as Bi{m'^^±,m'i.,,rn^,). 
The couplings are given in the following tabic: 
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Vk 


91 


92 


93 


Vl 


1 


W 


9H±i>ilL 


9w±w±v 


9w±hiL 


ul 


3 


w 


9H±ui^dL 


9w±w±v 


9w±ULdL 


ii 


bi 


TJ^± 


9H±iibi 


9w±w±v 


9w±hbi 


ii 






9H±iib2 


9w^w^v 


9w^hb2 




h 




9H±i2bi 


9w±w±v 


9w±i2bi 


i2 


h 




9H±Ub-, 


9w±w±v 


9w±Ub^ 



Counter-term and self energies diagrams: 

As explained in section. II-B, the counter-term A is fixed by real part of W^H^ self energy mixing eq^l Hereafter 



we list the fermion and sfermion 'Sf^r contribution to 'Shw ~ '^hw 



Diagram Fig.[21la: 



yS 



^Hwik") = ^{(mb{gw9H+9w9H)+mt{9w9H + 9w9H))Bi 



87r2 

mb{9w9H + 9w9H)BG)} 



(B19) 



Where k is the momentum of the external particle, gw,H = gwtb,Htb, the arguments of Bi functions are as 
Bi{P,ml,ml). 



F F 


9'h 


9§ 


9w 


9^ 


t b 


9Htb 
9HiO^± 


9Htb 
nR 


9wtb 

9wi('>^± 


9wtb 

9w^«x'^ 



Diagram Fig. |5]lb: 

^wC^^) = YQ:^9HS,Sj9wSiSj 
The arguments of i?i functions are as Bi{k'^ ^mg.^m^g ). 



Bo + 2Bi 



(B20) 



Si Sj 


9hs,s, 


9ws,Sj 


i>i 


II 


9h0iIl 


9 W 9111, 


ul 


dL 


9HuLdL 


9WuLdL 


ii 


bi 


9Hiibi 


9wiibi 


ii 




9Hii b2 


9wiib2 




bi 


9Hi2bi 


9wi2bi 


i2 


h 


9Hi2b,. 


9wh.b2 



Diagram Fig. Ellc: 



87r2 



((5i9G + 9§9h) Mm% ) + mp, i{gjig§ + g^g^) 



mp- 



{gnga + gHgG)TT^F, )Bo + {g^g'^ + g'haWB 



^L L 



L R 



(B21) 



The arguments of i?i functions are as -B(fc2, m|n. , to|t^. ). 





g'n 


9h 


9h 


9g 


t b 


9Htb 
9HyH± 


9Htb 
nR 


gctb 

9GiO^± 


n^ 
9Gtb 

9Gi0y± 



Diagram Fig. Elld: 



^HG — 



167r2 



gHS,Sj9GSiS, Bo{k'^, ml^ , ml. ) 



(B22) 
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Si Sj 


dHSiSj 


gCSiSj 


Vl 


II 




9g9iIl 


UL 


dL 


9HuLdL 




ii 


bi 


dHiibi 


dciibi 


ii 


h 


9Hiib2 


Set 162 




bi 




dclibi 


i2 


h 


9huU 


9Gl,b., 



The amplitudes of the counter-terms Fig. 2.2a — > d, are given by: 
Diagram Fig. |2]2a: 



= eA 



Diagram Fig.|l2b: 
Diagram Fig. |2l2c: 

Diagram Fig. |2l2d: 



-^2 — —gHHV^ 



gHGV 



m 



w 



T2 



gwwv 



gwwv 



-A 



(B23) 



(B24) 



(B25) 



(B26) 



where A is fixed by eg 1121 and the couphngs gnnvj 9hgv and gwwv have been defined above. 
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